The Electrostatic Field in Vacuum Tubes With Arbitrarily 
Spaced Elements 

By W. R. BENNETT and L. C. PETERSON 

VACUUM tubes with close spacing between electrodes have become of 
increasing importance in recent years. The higher transconductances 
and lower electron transit times thus obtained combine with other features to 
raise both the frequency and band width at which the tube may operate 
satisfactorily as an amplifier. Specific designs have been discussed in papers 
by E. D. McArthur and E. F. Peterson 1 , and by Fremlin, Hall and Shatford 2 . 
The important contributions to structural technique made by E. V. Neher 
have been described in the Radiation Laboratory Series 3 . Further im- 
portant advances in the art have been recently announced by J. A. Morton 
and R. M. Ryder of the Bell Laboratories at the recent I.R.E. Electronics 
Conference held at Cornell University in June, 1948. The material of the 
present paper represents work done by the authors over a decade ago, and 
naturally there has been considerable publication on related topics in the 
intervening years. It has been suggested by our colleagues, however, that 
some of the results are not available in the technical literature and are of 
sufficient utility to warrant a belated publication. These results have to do 
with the variation of the electric intensity, amplification factor, and current 
density which takes place along the cathode surface because of the nearby 
grid wires. 

We shall deal mainly with the approximate solution which neglects the 
effect of space charge. The correction required to take account of space 
charge is in general relatively small as shown by both qualitative argument 
and experimental data in an early paper by R. W. King' 5 . More recent 
theoretical work 19 extending into the high frequency realm has confirmed 
the minor nature of the modification needed. The problem is thereby re- 
duced to one of finding solutions of Laplace's equation which reduce to con- 
stant values on the cathode, grid, and anode surfaces. The original work 
on this problem was done by Maxwell 4 who calculated the electrostatic 
screening effect of a wire grating between conducting planes long before the 
vacuum tube was invented. All subsequent work has followed the methods 
outlined by Maxwell. In particular he suggested the replacement of the 
conducting planes by an infinite series of images of the grid wires and 
described an appropriate solution in series for the case of finite size wires. 
The useful approximation obtained when the diameter of the grid wires is 
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assumed small compared to their spacing was discussed in detail only for the 
case of large distances between the grating and each of the conducting planes. 
Figure 1 shows the assumed geometry of the grid, anode, and cathode. 
End effects are neglected. The origin is taken at the center of one of the 
grid wires which have radius c, and the X-axis is along the grid plane. The 
spacing of the wires between centers is a, the distance from grid to anode is 
</ 2 , and that from grid to cathode is d v No restrictions are placed on the 
sizes of a, do, and <k. Above the anode and below the cathode is shown a 
doubly infinite set of images which may be inserted to replace the conducting 
planes of the anode and cathode. By symmetry the potential from the 
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Fig. 1 — Array of images for production of equipotential surfaces in planar triode. 

array of charges there shown must be constant for all x when y = </« and 
also for all x when y = —d\. The double periodicity of the array suggests 
immediately an application of elliptic functions. The solution of the sym- 
metrical case was actually stated in terms of the elliptic function sn z by 
F. Noether 5 . The extension to the non-symmetrical case shown in Fig. 2 
is fairly obvious. One of the authors worked out such a solution in terms 
of Jacobi's Theta functions in 1935, but abandoned any plans for publishing 
his analysis in view of the excellent treatment appearing shortly after that 
time in the Proceedings of the Royal Society by Rosenhead and Daymond 6 , 
who applied Theta functions to both tetrodes and triodes, and both cylindri- 
cal and planar tube structures for the case of fine grid wires. Some of their 
formulas were later included in a book by Strutt 7 . Methods of calculating 
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the case of thick grid wires in terms of expansions in series of elliptic functions 
were discussed by Knight, Howland and McMullen 8-10 . The problem of a 
finite number of grid wires was treated by Barkas 11 . More recently tubes 
with close spacing between grid and cathode, but with anode and grid as- 
sumed far apart, have been analyzed in terms of elementary functions by 
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Fig. 2 — Variation of amplification factor along the cathode surface of a triode. 

Fremlin 12 . A solution based on the Schwartz-Christoffel transformation 
has been given by Heme 13 for the case of grid wires of finite size and approxi- 
mately circular in shape. 

Since the derivations have been adequately covered in the references cited, 
we merely state the final formula here and indicate how it may be verified 
as correct. Let V(x, y) represent the potential function corresponding to 
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Fig. 1, the planar triode with fine grid wires. The potential of the cathode 
is set equal to zero. Then in the space between anode and cathode, 

AV(x, y) - [2t d,(y - d 2 )/a + (d x + d*)f(x, y)]V, 

+ [B(y + di) - 2tt <ky/a - d x f{x, y)]V p , 



(1) 



where 



f(x, y) = (n 



t?i fr(* + *y - 2i rf 2 )/fl] 



(2) 



#i K* + iy)/a 
A = {dx + d,)B - 2* dl/a (3) 

atii(2iri di/a) 



B = In 



■kc$[{Q) 



(4) 



Here we have used Jacobi's notation for the d r function, as explained by 
Whittaker and Watson 14 , rather than the Tannery-Molk notation used by 
Rosenhead and Daymond. We write iMttz) for their &i(z). In our notation 

*G0 = 2 f: (-)V (B+1,a) * r sin (2» + 1)2 (5) 

n-0 

where the parameter t in the above formulas is given by: 

r = 2i(4 + 4)/« (6) 

By t?((z) is meant the derivative with respect to z: 

tfx'GO = 2 i; (-) n (2« + l) e ' (n+1/2),r cos (2» + 1) 2 (7) 

n-0 

Verification of the solution is straightforward. The resulting V(x, y) is 
seen to be the real part of a function which is analytic in the complex variable 
x + iy except for logarithmic singularities at the points where the Theta 
functions vanish. Hence V(x, y) satisfies Laplace's equation in two dimen- 
sions in the region excluding the singular points. Since the zeros of t?j(z) 
occur at z = mir + nirr, where m and n take on all positive and negative 
values as well as zero, the singular points of the solution are at 

x + iy = ma + 2in(di + d 2 ) — 2ida 

and ic + iy = ma + 2in(di + <£>) 

which coincide with the centers of the image circles of Fig. 1. The logarith- 
mic singularities represent line charges with the first set arising from a i?i- 
function in the numerator, yielding a positive charge, and the second set 
from the t? r function in the denominator giving a negative sign. The 
equipotential curves are approximately circular in the neighborhood of the 
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charges and hence V(x, y) gives a constant potential on the surface of each 
grid wire if the radius of the grid wire is small compared with the spacing. 
We may show by direct substitution that V(x, y) becomes equal to V p at 
all points of the anode and equal to zero at all points of the cathode. On 
the anode we have y = d 2 which, when substituted in the expression for 
f( x > y), gives the logarithm of the absolute value of the ratio of conjugate 
complex quantities, and hence 

fix, <h) = 

Substituting in (1), we then readily verify that V(x, d 2 ) = V p . On the 
cathode we make use of the quasi-periodicity of the t?i-function, as ex- 
pressed by 

*(*) = -« IOrr+to) iy 1 (z + Tr) > (9) 

to prove 

f(x,-d 1 ) = 2 ^, (10) 

from which it follows that V(x, — dj = 0. To show that all grid wires 
are at the same potential, we make use of the other periodicity of the 
i?i-function, 

*(S + 7T) = -<MZ), (11) 

which shows that 

fix ± ma, y) = f(x, y), m = 0, 1, 2, • • • (12) 

It remains to prove that V actually approaches the value V g in the neigh- 
borhood of the typical wire, which may be taken at the origin since the 
solution repeats periodically with the wire spacing. We let 

x + iy = ce (13) 

and assume c/a <K 1. Expanding in power series in c/a, we find that the 
first order terms are included in: 



f(c cos 6, c sin 6) = £n 



di(-2in di/ a) 



= B (14) 



The sign of the argument of the t?i-function in the numerator is of no conse- 
quence since it does not affect the absolute value. Substituting back in (1), 
we then find 

Lim V(c cos d, c sin 0) = V g (15) 

e/a—Q 

The solution is thus completely established. 
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The quantities in which we are specifically interested are electric field, 
amplification factor, and current density. The electric field is equal to the 
negative gradient of the potential function. The amplification factor is 
found by taking the ratio of partial derivatives of the electric field at the 
cathode with respect to grid and anode voltages. The current density may 
then be studied for any assumed operating values of grid and plate voltages. 

To calculate the gradient we note that since V(x, y) is the real part of an 
analytic function W(z) = V + iU, it follows from the Cauchy-Riemann 
equations, 

wto-W-i°V- -E x + iE y da) 

dx ay 

where E x and E v are the x- and y- components of the electric intensity. 
From (1), 

AW(z) = [(dx + <k) F(«) - 2*fr(« + <k)/o]V g 

+ [B(di - iz) + 2vi d 2 z/a - dx F(z)]V p 
where 

wo-fa *^rfr ^ {18) 

Calculating the derivative and making use of the relation, 

d[(z - xt) = diCO , 2i ( 19 \ 

di(z - ttt) &i(z) 

we find at the cathode surface 

F'(x-id l ) = — [1+C(*)1 (20) 

a 



17) 



where 



nl , T &iW(x + tdi)/a] / ?1 v 

C(x) = Im — . . — 7—^j\T~-\ ^ Zl) 

It follows that when y = - d x , we must have E x = and 

a^/2*- - [di + (A + d*}Cb)]V a 

+ [</ 2 - <*i - s aB/2v - diC{x)]V, 

The amplification factor is then given by 

= BEjSV g d x + (di + <fe)C(s) (23) 

M 3£,/3V, d 2 - di- aB/2ir - dxC{x) 

Numerical calculation from these formulas can be made by means of (5) 
and (6). When dx and (h are both large compared with unity, Eq. (23) 
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reduces to the familiar approximate formula derived, for example, in an 
early paper by R. W. King 15 , 

27Tfi?2 
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Fig. 3 — Variation of cathode field strength in a triode. 



Some calculated curves for /x and E„ are shown in Figs. 2 and 3. Figure 2 
shows the amplification factor as a function of the distance along the cathode 
with the ratio of grid-cathode separation to grid wire spacing as a parameter. 
The ratio of grid-anode separation to grid wire spacing is held constant at 
five. Only half the grid spacing interval is included since the curves are 
symmetrical. The increase in ^-variation as the grid-cathode separation 
becomes small is clearly demonstrated. For negligible /x- variation we must 
select d\/a of the order of 2 or greater. 
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Figure 3 shows the variation in field strength along the cathode for the 
typical operating point, V„ = -2 and V p = 100 volts. It is to be noted 
that for di/a less than 0.6, the electric field actually changes sign as we move 
from a point immediately below a grid wire to the midpoint between two 
grid wires. In other words a part of the cathode will not emit at all in these 
cases while the remainder emits in a non-uniform manner. In the rather 
extreme case of di/a = 0.4 only about a quarter of the cathode is emitting. 
It is worth noting how relatively rapid the "shadow" or "island" formation 
increases between di/a = 0.64 and 0.5 as compared to the increase in the 
interval from 0.5 to 0.4. 

If the equation for n is solved for C(x) and the result substituted back, in 
the expression for E v at the cathode we find: 

_ F . V + V P /n (24) 

v d x 4- (rfi + d,)/n 

where here of course /x varies with x. This is identical with the expression 
derived by Benham 16 from Maxwell's approximate solution except that in 
the latter case y. was a constant. Our colleague, Mr. L. R. Walker, has 
pointed out that the equation follows directly from the assumption of small 
grid wires without explicit solution for the potential function. Since the 
charge density <r on the cathode is proportional to the field strength (the 
factor of proportionality in MKS units is the dielectric constant e of vacuum 
or 9.854 X 10~ 12 farads/meter), Maxwell's capacity coefficients C gc and 
C pe may be calculated from 

a B = eE v = -(C 0C V„ + C pe V p ) (25) 

The minus sign is used here because we are taking the ratio of charge to 
voltage at the negative plate of the condenser consisting of cathode, grid and 
anode surfaces. Hence 

C " C = d X + (* + *)/M 

r - 6 A* * (27) 

Cpc "" <*i + (4 + *)/m 

Since p is variable, an integration is required to determine the total capaci- 
tance. From the periodicity of n with grid spacing it is possible to express 
the result in terms of the average values of C ac and C pc over an interval of 
length a along a direction parallel to the grid plane and multiply these 
values by the total area of cathode surface. 
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Equation (24) may be interpreted in a number of different ways of which 
we shall mention the following two: 

1. The "equivalent voltage" V + V p /n does not act at the grid but at a 
distance D from the cathode, where 

D = d x + (<*, + </,)//* (28) 

Both the equivalent voltage and distance vary along the cathode surface. 

2. The "equivalent voltage" 

V. = {V. + VM/[1 + (1 + <y«fr)/ M ] (29) 

acts in the grid plane and varies with distance along the cathode surface. 
As far as the cold tube is concerned the two formulas are equivalent at 
the cathode, but not at the grid. When the tube is heated and complete 
space charge is present, the two formulas also differ at the cathode. The 
current density in the presence of space charge is, according to (28) and 
Child's law: 



while, from (29), 



/ = K(V + VM*' 2 /D> (30) 



/ = KVy*/d\ (31) 



In both, K 2 = 32 e 2 e/81 m, where e/m is the ratio of electronic charge to mass. 
The value of current given by (31) is [1 + (1 + tk/dd/n] 112 times as large 
as that given by (30). If u > > 1 + d 2 /di the two values are nearly the 
same. In tubes with close grid-to-cathode spacing the inequality may not 
be fulfilled. As to which viewpoint is more accurate, we note that Ferris 
and North in their papers 17, 18 on input loading adopted the latter, and that 
at high frequencies where electron transit time must be considered the second 
viewpoint is preferable because of the more accurate representation of 
effects at the grid. For a more complete discussion see Reference 19. 
Figure 4 shows curves of relative current density as a function of distance 
along the cathode as computed from Eq. (31). The transconductance for 
unit area of cathode surface as computed from the same equation is given by: 

- 3.512(7, + Vp/nYIKdi/Dyi* micromhos. (32) 

The resulting variation with distance along the cathode is shown in Fig. 5. 
Defining the figure of merit M at a point x along the cathode as the ratio 
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between the transconductance dI/dV and the sum of C ge and C p0 at this 
point, we find from (30) 

M = (4//3) 1 "!*! + (rfi + < / 2 )/m]- 1/ V(m + 1) (M) 

where J = el /me, e/tn = 1.77 X 10 u coulombs/kg. From (31), we find 
on the other hand 

M = (47/3 diY V(M + 1) (34) 




O 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 
DISTANCE ALONG CATHODE ,*/a 

Fig. 4 — Variation of current density in a triode. 

Both formulas indicate that for a cathode capable of supplying a given cur- 
rent density the only means of improvement lies in decreasing the cathode- 
grid spacing. The improvement is extremely slow; doubling the figure of 
merit requires an eight-fold decrease in spacing. 

We again emphasize that the calculated current densities and figures of 
merit are functions of .v, the distance along the cathode. The total current 
between the two grid wires is found from (30) to be 

I T = 2K f 2 (V. + V p /n) 312 dx/D 2 (35) 
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Fig. 5 — Variation of Iransconductance along the cathode surface of a triode. 
while, from (31), 

It = £ f " [(mP. + Fj/Oi + 1 + h/dt)]*' dx (36) 

ai 2 Jx a 

where :v is given by 

V. + F p /m(.v ) = (37) 

On the basis of several reasonable assumptions it may be shown that both 
(35) and (36) lead to an approximate 5/2 power law instead of 3/2 power 
law. Such a law has actually been observed in cases where shadow forma- 
tion was suspected. 

We wish to express our appreciation to Messrs. R. K. Potter, J. A. Morton, 
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and R. M. Ryder for their encouragement, and to Miss M. C. Packer for 
aid in the numerical computations. 
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